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It is argued that, in nonextensive statistical mechanics with Tsallis en- 
tropy, the factorization of compound probability over subsystems is a conse- 
quence of the existence of thermodynamic equilibrium in the composite system 
and should be respected by all exact calculations concerning equilibrium sub- 
system. Using nonadditive energy satisfying this factorization, we propose 
an additive formalism of nonextensive statistical mechanics with additive q- 
deformed physical quantities and exponential distributions. This formalism 
leads to exact quantum gas distributions different from those given by fac- 
torization approximation with additive energy. The fermion distribution of 
present work shows similar characteristics to the distribution of strongly cor- 
related electrons given by numerical analysis with Kondo t-J model. 
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I. INTRODUCTION 

In this paper we will discuss some problems of nonextensive statistical mechanics (NSM) 
relevant to the factorization of compound probability into product of single body probability 

N 

P=Ilpn, (1) 

n=l 

where is the number of subsystems in the system of interest, p is the g-exponential dis- 
tribution (QED) : p oc [1 — (1 — q)/3H]~ ([.] > 0), as given by the maximization of Tsallis 
entropy Sq = — Tr^^ (Boltzmann constant A;^ = 1 and g > 0) under some constraints 



|T|-p. This factorization Eq.(|T|) has been viewed as a result of the independence of nonin- 
teracting subsystems having additive energy, just as in Boltzmann-Gibbs statistics (BGS) 
supposing short range interactions, and caused confusions in some theoretical studies of 
NSM and its applications to many-body problems. On the basis of a new idea relating 
Eq. (|T]) to thermodynamic equilibrium, we will argue that these confusions can be avoided if 
we introduce suitable nonadditive thermodynamic variables satisfying Eq.(|l|). Some theo- 
retical consequences for quantum distributions of this "equilibrium version" of NSM will be 
studied. 

Due to the necessity of defining additive average value of some extensive g-deformed 
thermodynamic variables, the discussions will be made within the formalism of incomplete 
statistics (IS) with Trp^ = 1 and normalized average x = Trp^x [Q. The reader will find 
that the quantum distributions of IS indeed show some particular properties already noticed 
with strongly correlated electrons. 

II. ABOUT INCOMPLETE NORMALIZATION 

IS as a alternative version of NSM was originally motivated by some theoretical pecu- 
liarities in the last Tsallis version of NSM based on the conventional normalization and 
unnormalized expectations 0. 

The basic assumption of IS is that our knowledge about physical systems is in general 
incomplete due to unknown space-time correlations or the effects of known interactions 
which can not be exactly studied (such as chaos). In this case, probability distributions are 
incomplete, i.e., Trp = Q 7^ 1 0] (or J2T=iPi = Q where w is only the number of accessible 
states in phase space). One can only write TrF(p) = 1 where F is certain function of p. In 
the case of complete or approximately complete distribution (such as in BGS), F is identity 
function. Recently, in order to overcome some of theoretical difficulties of NSM in keeping 
the framework prescribed by Tsallis entropy, we proposed F{p) = p'^ so that 

Trp^ = 1, (2) 
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where q is the incompleteness index 0. Since p < 1, we have to set g G [0, oo]. q = should 
be avoided because it leads to p = for all states. We note that Eq.(|^) has been successfully 
employed to deduce some power laws based on Renyi's entropy 0. 

This kind of incomplete normalization is possible whenever the phase space is partially 
known or accessible. With a fractal or chaotic phase space, e.g., a complete calculation of 
probability becomes in general impossible. In this sense, a plausible justification of Eq.@ 
may be inspired by a work of Tsallis ^ discussing nonadditive energy and probability 
distributions on fractal supports, although at that stage the work was not connected to 
anomalous normalization like Eq.(^. In that work, considering some simple self-similar 
fractal structures (e.g. Cantor set), one can obtain : 

where Vi{k) may be seen as the segments of the fractal structure at a given iteration of order 
k, V{0) a characteristic volume of the fractal structure embedded in a d-dimension Euclidean 
space, df = is the fractal dimension, n the number of segments replacing a segment of 
the precedent iteration, m the scale factor of the iterations and W = the total number 
of segments at the k^^ iteration. If we suppose that the fractal structure with k ^ oo 
is a phase space containing homogeneously distributed points, the exact microcanonical 
probability distribution of the k^^ iteration can be defined as 

^ V,{k) _ V,{k) 

where V is the total volume of the phase space. This distribution obviously sums to one. The 
problem is that V is an indefinite volume as k oo and impertinent for exact probability 
definition. In addition, V is not differentiable and contains inaccessible points. Thus exact 
summation in V would be impossible. Now if we define pi = ^^^^ as a physical or effective 
distribution, then we have 

w 

1=1 
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where Vq is a completely accessible and infinitely differentiable support on which the calcu- 
lation of Pi is possible. If we choose Vq = V{0), we can write Eq.(0) with q = df/d. The 
conventional normalization J2^iPi = 1 can be recovered when df = d. 

The above example is only a case of equiprobable distribution on simple fractal structure, 
but it illustrates very well the possibility that, in complex physical probability 

may not sum to one and may sum to unity only through a kind of power normalization, 
which, pertinent and useful for incomplete distributions, is consistent with the discussions of 
reference [|] on the mass calculation and the information consideration in porous structures. 

III. FACTORIZATION OF COMPOUND PROBABILITY AND 
THERMODYNAMIC EQUILIBRIUM 

In NSM, there are two major problems connected tightly to the factorization of compound 
probability. The first concerns the application of NSM to many-body systems via one- 
body distribution. NSM is originally intended to describe complex systems with long range 
interactions or fractal structure of space-time showing nonextensive phenomena. So from the 
beginning of this theory, Eq.(|l]) is supposed for composite systems containing statistically 
independent subsystems in order to elicit the nonextensive character by following relation : 

N 

ln[l + (1 - q)S,] = Y: ln[l + (1 - q)S,{n)] (4) 

ra=l 

For N = 2, Sg = Sg{l) + Sg{2) + {1 — q) Sq{l) Sq{2) as one often finds in the literature. Due to 
this independence, it has been believed by many that exact calculations within NSM should 
use the additive hamiltonian Hq = J2n=i ^n, where is the hamiltonian of n*^ subsystem 
Pf7|-[TT|. However, this hamiltonian is not compatible with neither Eq.(|l|) nor Eq.(P since 



these equations applied to QED mean ||1],[T2|,|T3[ : 

N N N k 

H=Y.H^+Y.[{q-i)pr' y: n^n. (5) 

n=l k=2 ni<n2<...<nfc j=l 

= Hq + He, 
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where j3 is the inverse temperature. In order to reconcile Hq and Eq.(|l|), a so called factor- 
ization approximation is proposed [|T^ by neglecting the second term at the right hand side 
of Eq.(^. This approximation has been, explicitly or not, employed in most of the applica- 
tions of NSM ITSHI^] via one-body QED. These applications certainly shows the usefulness 
of one-body QED, but the approximation neglecting the correlation energy by supposing 
sometimes weak interacting dilute particles [l^ is not a reassuring basis. Indeed, some re- 



cent works show that the correlation energy {He) given by the second term of Eq.(|^) is in 



general not negligible [0 and that the partition function given by using additive energy is 



completely different from that given by using Eq.(^ when is large |jTT|. So a doubt arises 



about the connection between independence of subsystems and additive energy. Recently, 
an interesting idea is forwarded to define a "quasi-independence" according to nonadditive 



energy Eq.(^ in order to apply NSM to turbulence flow problems |T^. As a matter of fact, 
this proposal implies rejection of classical independence for Eq.(|l]). 

The second problem connected to probability factorization is the establishment of zeroth 
law and the definition of temperature for NSM. It was believed that the zeroth law of 



thermodynamics was absent within NSM due to the paradox between Eq.(0) and the 



additive energy. Recently, a series of works have been published on this issue claiming 
the establishment of zeroth law and the definition of a generalized temperature on the basis 
of additive hamiltonian Hq and Eq.@ by neglecting H^. It is evident that the paradox 
mentioned above persists behind this approximate zeroth law. 

The central question is : Eq. (|l|) certainly implies independence of noninteracting systems 
for BGS, but does it mean the same thing for NSM? Very recently, Abe p2[ proposed a 



general pseudoadditivity for entropy required by the existence of thermal equilibrium in com- 
posite nonextensive systems. For a system containing A^ subsystems, the pseudoadditivity 
is : 

N 

ln[l + XsfiS)] = Y: ln[l + XsfiSn)], (6) 

n=l 

where / is certain differentiable function satisfying /(O) =0 and \s a constant depending 



on the nature of the system of interest. On the other hand, Eq.(P) apphed to TsaUis entropy 
means f{S) = S and A5 = 1 — g [Q, which directly leads to InTrp^ = Y^iLi InTrp^ or Eq.(|I]) 
(i.e. with classical probability pi of the state i, {piPjY = pjj means piPj = pij). So Eq.(|l]) has 
nothing to do with statistical independence of subsystems. It is a consequence of the existence 
of thermodynamic equilibrium and must be rigorously respected by all exact calculations. 
Equilibrium energy has been proved ||23| to satisfy the same kind of pseudoadditivity as 



Eq.® ("S" is replaced by H). If we choose f{H) = H and A// = ((? — we get 

N 

ln[l + (g - 1)PH] = Y: ln[l + (g - l)PH^] (7) 

n=l 

which is just Eq.(|^) satisfying Eq.([l|). In this way, the zeroth law becomes evident and a 
temperature can be straightforwardly defined at maximum entropy and minimum energy 

10- 



IV. ADDITIVE FORMALISM OF NSM 



A. Information measure 



The (/-logarithmic information measure 



- (8) 

is a nonadditive generalization of Hartley formula / = ln(l/p) and can be employed to 
deduce TsaUis entropy |]I],|3|,P^. Ig or / is the information needed to specify at which state 
the system is localized, v equals 1 — g or g — 1, depending on the normalization procedures 
of p [^. Using Eq.(p, we get : 

N 

ln(l + z//.) = ^ln(l + z/4")) (9) 

n=l 

where J^"^ is the information needed to specify the n*^ subsystem. This pseudoadditivity is 
evident if we recast the generalized Hartley formula Eq.(D as follows 

/. = = (10) 

V V 



where 

,.i„l = Ml±i^, (11) 

P ^ 

can be referred to as g-deformed information measure and is additive supposed Eq.(|I]). It 
is noteworthy that this I is not the quantity of Hartley information if p is a nonextensive 
distribution for v ^ {]. 



B. Canonical ensemble 



Now let us define an additive entropy 5* as follows : 



P 

and an additive g-deformed "hamiltonian" 



S = lTpnn- (12) 



_ ln[l + (g - l)PH] 



So Eq.(l^) becomes 



N 



h=Y.hn. (14) 



n=l 



This means following transformations : 



(g-l)/3/i _ 1 (g-l)/3fe„ _ 1 



and 



p = i[l + (g - = }-e-P\ (16) 

where Z"? = Tre~'^'* 0. It should be noticed that, when addressing a system of particles, 
we have to write if„ = |^ + \4 for single particle so that hn = '"^^"'"^'^ (g-^)^'""^^"^^ where ^ is 
the classical kinetic energy and Vn is the potential energy. It is clear that Hn, instead of hn, 
is the physical energy. When g = 1 {He = 0), we recover Hn = hn and H = h = J2n=i im- 
The g-deformed internal energy u is defined as follows : 



u = Trp'/i. (17) 



We can easily show that the distribution Eg . ([16|) can be yielded by the maximum of the 
additive "entropy" 5* (which surely exists due to the monotonic relation between / and 
lu) under the constraints of Eg. (pTf) and incomplete normalization Trp'^ = 1. It is easy 
to verify that S = InZ + jSu and, via the zeroth law, |^ = /? = 1/T. The "first law" is 
given by du = TdS — pdV where p is g-deformed pressure and V the volume of the system 
which is chosen to be additive here. The g-deformed Helmholtz free energy / is defined as 
f = u — TS = —T In Z and can be connected to the nonadditive one Fg = —T ^\^^^ [|1],|^ 
as follows : 



hl|l + (q- 1)I3F,] 

f fe^T)^ ■ '1*' 

So p = — (ff )t = P/Z^^'^ where P = — (§y-)r is the real pressure. In this scenario, 
the thermodynamic eguilibrium of a system C containing two eguilibrium systems A and B 
satisfying V{C) = V{A) + V{B) corresponds to P{A) = [3{B) andp(A) = p{B). This implies 
that P{A) 7^ P{B) if Z{A) ^ Z{B). This is because we have supposed nonadditive energy 
and additive volume. As a matter of fact, in this formalism, if we want that P{A) = P{B) 
at eguilibrium, we must accept nonadditive volume and additive g-deformed volume v with 
which the first law becomes du = TdS — Pdv. This means : P = — (^)t = ~(|f )t- This 
relation can help to deduce v — V relation. We also have v = {§p)t where the g-deformed 
Gibbs energy is given hj g = f + Pv. Since v is additive, V will be nonadditive if v is not 
a linear function of V. We will come back to this issue later in this paper. 



C. Grand canonical ensemble 



It is known that the grand canonical ensemble QED has been given by 

poc[l-(l-g)/?(//-/xiV)]T^ (19) 

for N identical particle systems, where fi is chemical potential. This distribution has been 
widely used for guantum particle systems |l^- |l6| , ^ . But the zeroth law has never been 
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rigorously established for this ensemble. In a previous work ||26|, one of the authors of 
present paper deduced exact quantum distributions on the basis of Eq.(|19|) and following 
relation suggested by Eq.(|T]) : 

N 

pcx[l-(l-g)/?(if„-/i)] — . (20) 
In the framework of IS [0,^], the exact distributions are given by 

nu = s (21) 

[l + (g-l)/?(e,-/x)]-±l 

where is the energy of one-particle state k and"+" and "-" correspond to fermions and 
bosons, respectively. Now we show that this distribution can be written in exponential 
form just as for conventional noninteracting quantum gases and that the zeroth law can be 
rigorously verified. 

Let us suppose /? = = - (g - !)/?>'] (or (3' = and ^i = 

(or fi' = = ~ {q — 1)/?/^]) which imply P'fi' = Eq.(|2D|) can be recast as : 

1 N N 

p = - [1 - (1 - q)P'Kn] — [l + (1 - g)/9V]- (22) 

- ]_p-NI3{h„-u;) 

z 

where co = "^'^l/l^gf and = mfi+fa-iyH^] ^ ^ ^^^^ _ _ ^^^^^ _ ^^^^^ ^ 
{Tr[l - (1 - q)f3'Hn]^[l + (1 - g)/5y]^}^ = (Tre-«^(^"-^))^ = and where z is one- 
particle partition function. Just as for canonical ensemble, this exponential distribution can 
be shown to be the result of the maximization of S under the constraint N = Trp'^N in 
addition to Eq.([T7|) and incomplete normalization. Now Eq. (pT]) can be written as 

1 1 

Uk = Trp^Uk = — T"— = q 5 (23) 

where is the eigenvalue of hn- From Eq. (|2TD , we see that, for free particles {in the sense 
that we do not write the energy of interaction between particles in the hamiltonian and let it 
he "absorbed" in the nonextensive part of energy He and related to q different from unity), 
we have to set g < 1 to ensure positive [1 + (g — l)(3{Hn — p)] for fermions when T ^ 0. 



This means that, at low temperatures, there will be few fermions beyond Fermi-energy. For 
bosons with yU < 0, we have to put g > 1. 

It is straightforward to write S = InZ + (3'u + (3'ujN and, with the method employed in 
references [^|12|) to show that, for a system C containing two equilibrium systems A and 
B satisfying N{C) = N{A) + N{B), (3'{A) = (3'{B) and uj{A) = uj{B) thus ij'{A) = ^\B), 
^{A) = fi{B) and P{A) = f3{B) or T{A) = T{B). The zeroth law is verified. One may ask 
why we identify /3, instead of to real temperature. The possible reasons are the following 
: 1) (3 is the Lagrange multiplier of the constraint on real energy in entropy maximization; 2) 
if (3' is inverse temperature, then /i' must be chemical potential, which makes it impossible 
to get distribution Eq.(^) by real entropy maximization with the constraint on iV; 3) e/ 
would be different from the chemical potential fx' and equal to ij^^i'^q^p, ^, which inevitably 
drops to zero when T ^ 0. 

Now we will focus the discussion on 2D fermion distribution. According to Eq.(|21|) the 
Fermi energy e/(, at T = is given by e/^ = where a is the particle density and h 

Planck constant. 

When T > 0, the summation N = J2k^k can not be calculated for arbitrary q to give 
explicit expression of ej. So we have recourse to numerical calculation of g-dependence of 
for given temperatures (Figure 1) and T-dependence for given q values (Figure 2). We see 
that, for the approximate distributions functions (ADF) n/c = l/{[l + (q'— l)/5(efc— + 
deduced from Eq.([l^) with factorization approximation and additive energy |T^, ej depends 
only slightly on q. On the other hand, the ej of IS in the present work shows a strong increase 
with decreasing q up to two times ejg of conventional Fermi-Dirac distribution (CFD). This 
Cf increase has been indeed noticed through numerical calculations for strongly correlated 
heavy electrons on the basis of tight-binding Kondo lattice model p7|j28|] . In Figure 2, we 



show that of IS does not monotonically decrease with increasing temperature, in contrast 
with the 6/ behavior of CFD. This kind of non monotonic temperature evolution of e/ was 
also seen through numerical work for correlated electrons in two-dimensional t — J model 
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The IS distribution given by Eq.(^) is plotted in Figure 3 for T = lOO-ft'. The particle 
density o is chosen to given c/q = 1 eV. As expected, the distribution changes drastically 
with q. When q decreases, we notice a flattening of the distribution with always a sharp drop 
of occupation number at e/ which increases. This flattening is also noticed in numerical 
calculations for strong coupling electrons E^. 



The grand canonical partition function Z can be calculated to give the g-deformed grand 
potential uj = —T' InZ = ^ J2k 1^(1 ~ ^k) as usual, where T' = = T + (1 — The 
grand potential is given by : 



(,-1) ■ (^^' 

In this g-deformed extensive formalism, Euler theorem applies just as in BGS. So we 
have g = ujN = u + pV — T'S. Compare this with S = InZ + (3'u — (3'ujN to obtain the 
following equations of state : 

T' 

pV = T'\nZ = —Y,^n{l-nk). (25) 

^ k 



V. DISCUSSION AND CONCLUSION 

The formalism of NSM presented here is required by the existence of thermodynamic 
equilibrium or by Eq.(|lD for nonextensive systems described by Tsallis entropy. Theoretically 
the formalism is self-consistent. Experimental or numerical evidences are needed to verify 
the thermodynamic relations. In this framework, all the successful applications of NSM 
conforming with Eq.(|I]) are still valid. But the approximate applications carried out for 
many-body systems using additive energy as exact hamiltonian (not consistent with the 
existence of thermodynamic equilibrium) should be carefully reviewed. 

As mentioned above, we have noticed similar properties between the IS fermion distri- 
bution and that of strongly correlated electrons |2^J2^. This similarity shows the merit of 
NSM in describing strong interacting systems. On the other hand, we noticed that a flat- 
tening of n drop at e/ with increasing correlation, observed experimentally and numerically 

11 



with weakly correlated electrons p7|,p8|,|30|-p^ , is absent within NSM fermion distributions 



which always show very sharp n drop at e j at low temperatures. A detailed study of this 
problem will be presented in another paper of ours. 

It is worth mentioning that, in the present work, energy is nonadditive to satisfy the 
requirement of thermodynamic equilibrium. Nonadditive energy can happen if interaction is 
no longer of short-range and not localized only between the containing walls of subsystems. 
But in the literature, there are rarely explicit expressions of nonadditive energy. One of 
explicit examples is the long-range Ising model ||33| where the internal energy U {N, T) may 



proportional to A^'^ (c is a constant), instead of A^, the number of spins in the system. This 
energy can be shown to satisfy Abe's pseudoadditivity for energy p3| if we choose f(U) to 
be proportional, e.g., to A^ (with Xh = 0) or to — 1 (with \h = 1). 

Indeed, theoretically, nonadditive physical quantities (energy, volume ... ) is not evident 
within the statistics with complete distributions because all possible states (all points in 
phase space) are counted and summed here. But from the viewpoint of incomplete statistics, 
nonadditivity may be interpreted as a consequence of incomplete summation of state points 
in phase space due to the incompleteness of our knowledge about the physical systems 0. 

Nevertheless, the fact that the correlation energy of NSM depends on temperature, as 
shown in Eq.(^ or Eq.(^, is not an easy aspect to be understood. A possible interpretation 
is that these nonadditive equations are required or prescribed by thermal equilibrium with 
Tsallis entropy and so naturally change with temperature. This implies that the effect of 
correlations may depend on temperature. 

Summing up, within the framework of incomplete statistics, it is argued that the nonex- 
tensive thermostatistics should be based on the factorization of compound probability sug- 
gested, not by "independence" of noninteracting systems, but by existence of thermody- 
namic equilibrium in interacting systems having Tsallis entropy. So this factorization must 
be viewed as a fundamental hypothesis of NSM and rigorously satisfied by all exact cal- 
culations relative to equilibrium systems. On this basis, we have elaborated an additive 
formalism of NSM based on the maximization of an additive deformed entropy subject to 

12 



constraints on additive particle number and additive g-deformed energy. The IS quantum 
distributions of this formahsm are compared with the distributions previously obtained by 
using complete probability and additive energy in factorization approximation and also with 
numerical results for strongly correlated electrons. It is shown that some effects of strong 
correlations : the flattening of the fermion distribution and the sharp cutoff of occupation 
number at e/ which shows strong increase with increasing interaction, can be observed in IS 
fermion distribution with decreasing q value. 
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Figures: 

FIG. 1. g-dependence of Fermi energy e/ of IS quantum distribution in present work and of ADF 
given by the factorisation approximation with additive energy. The fermion density a is chosen to 
give e° = 1 eV for CFD distribution at T = 0. IS e/ shows strong increase with decreasing q up 
to two times ef^. But ADF e/ depends only slightly on q. We also notice that the T-dependence 
of IS e/ is not monotonic as shown in Figure 2. 

FIG. 2. T-dependence of Fermi energy of IS quantum distribution in present work. The 
fermion density a is chosen to give e° = 1 eV for CFD distribution at T = 0. The T-dependence 
of e/ is in general not monotonic, in contrast with the classical decreasing behavior of e/ with 
increasing temperature. We notice that, at low temperature, shows an increase with increasing 
T. 

FIG. 3. Fermion distributions of ADF and of IS in present work. ADF distribution is only 
slightly different from that at q = 1 (CFD one) even with q very different from unity. The IS 
distribution of present work changes drastically with decreasing q. As q ^ 0, the occupation 
number tends to 1/2 for all states below e/ which increases up to 2 times e/g. 
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